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We introduce higher dimensions into the problem of Bose-Einstein condensates in a double-well 
potential, taking into account orbital angular momentum. We completely characterize the eigen- 
states of this system, delineating new regimes via both analytical high-order perturbation theory 
and numerical exact diagonalization. Among these regimes are mixed Josephson- and Fock-like 
behavior, crossings in both excited and ground states, and shadows of macroscopic superposition 
states. 
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I. INTRODUCTION 

Bose-Einstein condensates (BECs) in double-well po- 
tentials continue to receive much attention due to their 
wide variety of applications, ranging from precision mea- 
surements [ij-dl to optical information processing [|| and 
quantum computing Being a natural realization 

of a macroscopic two-state system, they provide an ideal 
system for studying fundamental quantum many-body 
phenomena and Josephson- type effects [1, Q. BECs in 
double wells exhibit distinct physical regimes. Three pa- 
rameters are frequently used to distinguish them: the 
number of atoms N, the tunneling coefficient J, and 
the interaction coefficient U. The two main regimes 
are called the Josephson regime and the Fock regime. 
In the former, tunneling dominates over interactions, 
C/N > 1 with £ = J/\U\, and the limiting case is the non- 
interacting gas, f/ = 0. In contrast, in the Fock regime 
( < 1, interactions are much bigger than tunneling, and 
the limiting case is the infinite-barrier or zero-tunneling 
case, J = @,@|. These regime designations are based 
on just two single particle states, one in each well, and 
most approaches in the literature are an extension of this 
two- mode concept. We will show that when more single- 
particle states are taken into account many new regimes 
occur. Moreover, the dimensionality of the double- well 
potential manifests in the form of a new quantum num- 
ber, the orbital angular momentum. 

In this Article, we use a two-level generalization of 
the Lipkin-Meshkov-Glick (LMG) Hamiltonian 0, 03 
to investigate the behavior of ultracold bosons in a three- 
dimensional (3D) double-well potential. In a previous 
study, we relaxed two assumptions commonly made in 
similar systems: the s ymm etric-trap assumption and the 
one-level assumption [12l. fl3j. We showed that, for rea- 
sonable physical parameters, one already requires the 
first excited state on each side of the double well for the 
order of 10 to 100 atoms in typical BEC experimental 
systems. Therefore, in this study, we consider the ef- 
fects of this excited state in detail. The two levels of our 
3D Hamiltonian give rise to eight single-particle modes, 
and atoms in the upper or excited level are allowed to 
have nonzero orbital angular momentum, as illustrated 




FIG. 1: (Color online) Schematic of the potential, (a) Equipo- 
tential surfaces of the double-well potential and (b) zoom of 
single well and single particle eigenfunctions, where I is the 
orbital angular momentum and m its z projection. The first 
excited level is triply degenerate in 3D. 



in Fig. [T] We give the criteria that permit one to distin- 
guish between Josephson and Fock regimes in both levels, 
i.e., we consider also the hopping parameter and inter- 
action parameter for the excited level, leading to novel 
regimes. We characterize the eigenvectors and eigenval- 
ues in all regimes. In some regimes different kinds of 
macroscopic superposition (MS) states, i.e. states for 
which atoms simultaneosuly occupy both wells, are en- 
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countered, ranging from MS states involving atoms only 
in the bottom level to MS states between atoms only in 
the excited level with angular momentum, and including 
mixed MS states with atoms in both levels. We will illus- 
trate these different kinds of states with surface plots of 
their probability amplitudes as a function of both energy 
and their Fock index. The Fock index orders states in 
Fock space, as shall be described; the index increases as 
more and more atoms occupy the upper level. Pertur- 
bation theory can give rise to shadows which replicate 
unperturbed patterns in these plots, and appear as faint 
copies at higher Fock index. We give the criterion to der- 
macate the regime in which these shadows of MS states 
occur. 

Ultracold bosons in double wells are a good candidate 
for the experimental realization of MS states, and there 
are many theoretical proposals in this direction fl4l - |2"lj ]. 
One of the main reasons the MS problem has been so 
heavily pursued, besides technological applications [H-Q]; 
is because of the possibility of a breakdown in the predic- 
tions of quantum mechanics at a macroscopic level Q. 

In our investigations, we utilize two main methods: 
numerical exact diagonalization and analytical perturba- 
tion theory. We clearly delineate the different regimes 
in our two-level eight-mode Hamiltonian modeling a 3D 
double well. Our approach can be extended in a straight- 
forward manner to one- or two-dimensional systems with 
four or six modes, respectively. To illustrate the regimes 
we give numerical examples of every different case us- 
ing an experimentally realistic double well potential, the 
Duffing potential. We focus on statics, leaving dynami- 
cal considerations, such as tunneling in different regimes, 
for future work. Our presentation is organized as fol- 
lows. In Sec. fll] we describe our model and elucidate 
the energy scales relevant for the problem. Once the pa- 
rameters that determine the different regimes are clearly 
stated, we characterize the eigenvectors for the different 
regimes in Sec. IIIII In Sec. IIV1 we find criteria for the 
boundaries between different regimes and for the valid- 
ity of the one- and two-level approximation. In Sec. [V] 
we illustrate numerically the theoretical results and find 
expressions showing how all criteria vary as a function 
of the number of atoms N. In Sec. IVII we summarize 
and discuss our work. We relegate a detailed description 
of the construction of single-particle eigenstates and our 
perturbative methods to appendices. 



II. QUANTUM TWO-LEVEL APPROXIMATION 

A. Qualitative discussion of physical regimes 

The one-level or two-mode assumption utilized in the 
original LMG model is valid if coupling to higher single- 
particle energy levels can be neglected. Although levels 
are only completely decoupled when there are no atom- 
atom interactions, coupling to excited states is very small 
when the interactions are much weaker than the single- 



particle energy-level spacing for each well, AE. However, 
even in this regime effects of the excited level are still im- 
portant. Eigenvalue crossings, which cannot be described 
by a one-level approximation, occur when either the num- 
ber of atoms or the interaction energy is greater than a 
critical value [l2l . [l3j . There are two key crossing regimes 
in the Fock regime. First, states with definite occupation 
of the excited level emerge among the lowest-lying N + l 
eigenstates. Second, as the interactions are increased fur- 
ther, such crossings involve not only excited states but 
also the ground state; i.e., when the system is in the 
ground state one has a finite probability of measuring an 
atom in an excited level. In Sec. IIV1 we quantitatively 
delineate these regimes, and we show that the former 
crossing does not occur if the condition iV 2 |£/| <C A£J is 
fulfilled while the latter does not occur if N\U\ <C AE. 
These two key crossing regimes occur also in the Joseph- 
son regime for very small interactions, or even vanishing 
interactions, as detailed below. Finally, we use different 
criteria to clearly identify the regime of parameters for 
which the model is valid. In particular, our approach 
is valid under the dilutcness condition, or insofar as the 
interactions are not too strong. 

The very different physical phenomena found in dif- 
ferent regimes justify the range of theoretical approaches 
encountered in the extensive literature on the double- 
well problem. Joscphson oscillations were first predicted 
for the limiting case of a non-interacting g as, called 
the "extreme Josephson" or "Rabi" regime [22|]. For 
non-zero interaction but still in the Josephson regime, 
mean field (MF) app roaches have successfully predicted 
Josephson effects [23M25I ] , while as the atom interaction is 
increased these models predict macrosco pic self-t rappin g 
of the condensate in one of the wells [23l l24l . [26l - l28| . 
While a number of these works also present the quan- 
tum analysis of the system, a complete quantum-phase 
picture of the problem showing the correspondence with 
the MF approaches, and the transition from delocalized 
to a fully quantum regime was developed by Mahmud 
et al [2!|. Hence, the MF approach finds its limitation 
in the Fock regime, and two-mode approaches find their 
limitations insofar as more highly excited states are re- 
quired to describe the dynamics of the problem. The 
latter difficulty can be overcome using a multi-mode ap- 
proach, and thus extending the results to larger interac- 
tions [27j . but in general, for both cases other methods 
are required. One such method is multi-configurational 
time-dependent Hartree (MCTDH) theory [3(|, and its 
stationary counterpart [3l|, HII- Other approaches are 
based on exact diagonalization of the LMG Hamilto- 
lll. with more than two modes and additional 



12L 13l [33| | . The tunneling dynamics of the sys- 
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tern have been studied using MCTDH methods, finding 
different dynamics than MF methods [HI, [HJ]. In ad- 
dition, LMG methods predict exponentially long tun- 
neling times whereas MF methods predict macroscopic 
self-trapping poj . Tunneling and macro scop ic self trap- 
ping have been observed in experiments |36t |37| , though 
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the latter phenomenon can be attributed to long tun- 
neling times I38|. |39i . small asymmetries in the double- 
well potential [20 ] . or even inhomogeneities in the inter- 
actions fiof . 

Furthermore, the transition from the Josephson regime 
to the Fock regime, or from a coherent to an incoherent 
regime, shows the limitations of the MF approach [4lj . 
The study of this transition elucidated that, in the lat- 
ter regime, strongly-correlated quantum states appear, 
showing macroscopic occupation of two single particle 
states localized at each well [H|,[2!|[42j]. These states take 
the form of MS states, colloquially called Schrodinger cat 
states, and require two independent macroscopic modes 
that can be coupled or entangled. MS states were also 
proposed for two-species BECs, i.e., for two internal 
modes [43l - l47j . A two-s pec ies BEC in a single well is 
mathematically identical [48[ to a single-species BEC in a 
double- well in the one-level approximation. However, our 
use of angular momentum makes the double well problem 
quite different from the usual two-species BEC one. 

As atomic interactions arc increased, not only MS 
states appear, but also coupling to states with occupa- 
tion of the excited levels should be considered, this being 
responsible for fragmentation of the condensate (29l l49j . 
Therefore, in this regime, more levels should be included 
in the model; MCTDH provides a superior method, but 
is computationally limited compared to our approach; 
specifically, to-date MCTDH has not been able to treat 
3D systems with orbital angular momentum. 

To make an analogy, helpful for the general reader who 
may be familiar with the language of ultracold bosons 
in optical lattices [50l - l52| . suppose we considered not 
two wells but an infinite number of wells along a line. 
This is the lattice problem, and the LMG model be- 
comes the Bose-Hubbard model. Then our upper level 
orbital becomes the second band, which is D-fold degen- 
erate in D-dimcnsions. Fragmented states become the 
Mott-insulator, and the Fock regime is identified with 
the Mott-insulating regime, while the Josephson regime 
is identified with the supcrfluid regime. However, the 
Bose-Hubbard model is generally treated for small fill- 
ing factors of a few atoms per site, while our double- well 
problem is geared towards large filling factors in hopes of 
achieving an MS state. 



B. General Hamiltonian and double well potential 

The second-quantized Hamiltonian for a system of N 
interacting bosons of mass M confined by an external 
potential V(r) at zero temperature is given by 

h 2 



H= Jd 3 r&(r) 
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on the s-wave scattering length a s of the atoms, g = 
Airh 2 a s /M. 

We consider a 3D double- well potential with minima at 
r = ±r m ; n € K 3 and a local maximum at r = 0. Without 
loss of generality, we will consider a separable potential 
V(r) = V x (x) + V y (y) + V z (z), built up from two harmonic 
single-well potentials V x (x) + V y (y) = \{u 2 x 2 + uj 2 y 2 ), 
and a generic ID double-well potential in the third coor- 
dinate z, with two minima at z = ±z m i n /2 and a max- 
imum at z = 0. From here on we call the difference 
between the maximum and the minima in such a poten- 
tial the barrier height, denoted by Vo- Near a minimum, 
this ID potential is V(z ± z m m) ~ ^^z 2 where u> is an 
effective single-well trapping frequency. 

Equation (TIJ is valid at low densities, when only bi- 
nary collisions are relevant, and at low energies, when 
these collisions are characterized by the s-wave scatter- 
ing length of the atoms [53|. The diluteness condition 
for a weakly interacting Bose gas is y/\n of] -C 1, where 
n is the average density of the gas. In the context of 
the double-well potential, an upper bound on the den- 
sity of the gas is approximately n = N/(^/2tt aho) 3 , where 
dho = y/H/Mui is the oscillator length and w is the single- 
well harmonic oscillator frequency. Correspondingly, we 
restrict our discussion to the regime 
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Although the system is said to be "weakly interacting" 
when condition ^) is met, the interaction energy can 
be on the order of the kinetic energy, or even bigger as 
corresponds to the Fock regime @, Q- In Sec. IIVI we 
obtain condition ^ in terms of the relevant coefficients of 
the double well problem, thus permitting us to compare 
this criterion with the criterion characterizing the Fock 
regime in the numerical results given in Sec. [Vj 



C. Two-Level approximation 

Double-well potentials in one and two spatial dimen- 
sions can be achieved in extremely anisotropic traps. The 
ID and 2D transverse trapping frequencies must be suf- 
ficiently high to reduce the dimensionality of the single- 
particle wavef unctions, but should not be near any po- 
tential resonances 



54[. In this Article, we restrict our 



where ^(r) and ^(r) are the bosonic annihilation and 
creation field operators. The coupling constant g depends 



attention to the 3D case, in particular the axially sym- 
metric one, uj x — uj y . 

We can expand the field operators in any basis of the 
Hilbert space. We use a fixed single-particle basis, con- 
structed from the delocalized eigenfunctions of the sin- 
gle particle Hamiltonian H sp = — |gV 2 + V(r). Our 
site-localized basis is constructed from appropriate su- 
perpositions of delocalized eigenfunctions, analogous to 
how Wannier states are obtained from Bloch functions 
on a lattice [55|. This approach results in spatial states 
of form VWm(r — ij-), where j signifies the left or right 
well, n is the single-particle energy level, £ is the orbital 
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angular momentum in 3D, and m is its z-projection, as 
sketched in Fig. [T] see also App. A and Fig. |H] for a more 
detailed description. Then the field operators can be ex- 
panded in this basis as 



(3) 



where ri = — r min and r 2 = r m j n are the minima of the 
left and right wells. The operators &]„^ m and b^ nlm sat- 
isfy the usual bosonic annihilation and creation commu- 
tation relations, 
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£m > j n £ m 



0. 



(4) 



For aho <C z m in the functions VWm(r) closely resemble 
the eigenfunctions of the harmonic oscillator potential 



V(r) = ±w 2 r 2 : 



VWm(r) ~ Rnt{r)Y im {6 , </>), 



(5) 



for n £ {0, 1,2,.. .}, £ e {n, n — 2, n — 4, . . . , £ m i n }, and 
m e {-£, -l + l, . . . , t-l, £}. Here R n i{r) is the 
radial part of the wavefunction, Yi m (9, <p) are the familiar 
spherical harmonics, and £ m ; n is when n is even and 1 
when n is odd (56|. The energy of an atom associated 
with the wavefunction ip n £ m (r — ry) is 
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hio(n + 3/2). 



(0) 



We emphasize that the harmonic-oscillator description 
is only approximate; actual eigenfunctions are distorted 
from spherical harmonics as sketched in Fig. [TJ The two- 
level approximation, i.e., truncating n at 1, is the lowest 
order of n at which the dimensionality of the double- 
well becomes apparent. Because n = £ for n £ {0, 1}, 
both the total orbital angular momentum of an atom and 
its energy level are described by the quantum number 
£. In the two level approximation, the subscript n is 
superfluous and is hereafter suppressed. 



D. Two-Level Hamiltonian 

Substituting Eq. ([3]) into the second-quantized Hamil- 
tonian (fT]) yields the two-level Hamiltonian 



H — H + H i + Hqi , 



(7a) 



where 
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for j € {1, 2}, £ € {0, 1}, m,m' e {-£, £} and 
m" € { — 1,0,1}. Here the number operator fiji m = 

bjtmPjim- The one-level Hamiltonians H and Hi de- 
scribe atoms in the lowest and first excited energy levels, 
respectively Atoms in different energy levels are coupled 
by the operator H i. I n addition to the level spacing 
AE = Ei — E , the problem is characterized by the fol- 
lowing energies, 
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Three basic processes characterize the two-level Hamil- 
tonian: single-atom tunneling between left and right 
wells, two-atom hopping between energy levels, and 
atom-atom interactions. Individual atoms tunnel be- 
tween wells with energy Jg m . Such transitions do not 
alter the z-component of an atom's angular momen- 
tum m as we have chosen w x = uj y . Furthermore, 
single-atom transitions between energy levels £ are for- 
bidden by the orthogonality of the localized wavefunc- 
tions "i/^ m (r ± r m i n ); that is, we have chosen a basis such 
that only interactions can mix levels. Since one-atom 
hopping is related to the first term in Eq. ([T]), after ex- 
pressing the field operators in terms of the localized single 
particle eigenfunctions, all integrals involving pairs of lo- 
calized functions with different values of £ and m vanish. 
Instead, the second term in Eq. (flj gives same-site inter- 
level hopping, which is achieved by pairs of atoms. No- 
tice that every integral of the form J \ipj£m\ 2 ipje'm.''<Pj£"m" 
vanishes, as can be easily shown performing each inte- 
gral on variable 9 using the expressions for the local- 
ized functions detailed in App. [A] The only exception 
is for £' = £" = 1 and ml = 1, m" = —1, due to the 
fact that ipj£- m = ipjim- These give the hopping terms 
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h.c. 



and 
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that appear in Eq. ((TJ. Then, the only terms that have 
been neglected to obtain Eq. (TTJ correspond to off-site in- 
teractions. These correspond to terms whose integrand is 
of the form \ijjjt m \ 2 \ipj'£ m \ 2 or ijjj im ipj>£ m , which are simi- 
lar to those considered in [4^] . These integrals bring in a 
term proportional to exp [— (z m in/aho) 2 ] . Then, the cor- 
responding interaction coefficients are much smaller than 
any interaction between atoms in the same well 
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which do not show this term, if the barrier height, the 
distance between wells, or both are not too small. In par- 
ticular, the barrier height is big enough to neglect these 
terms in the Fock-like regimes that capture our interest 
for MS states. 

Thus, according to Eq. ((J), level hopping is induced by 
atom-atom interactions in the same well, unlike tunnel- 
ing between wells. Two interacting atoms hop together 
between energy levels with energy U^„ . A pair of atoms 
can hop from the lowest to the first excited energy level 
of the same well in two distinct ways. Either both atoms 



enter the 



= state of the excited level or one atom 



enters the to" = +1 state and the other the m" = — 1 
state. A similar process, described by the interaction en- 
ergy Uli, occurs within the excited level. Pairs of atoms 
in the same well interact with energy Up 1 ^,. 

Here we have used eight modes to expand the field 
operators in three dimensions. The formalism up to this 
point is equally valid if we consider more levels, i.e., more 
modes, yielding to new terms in the Hamiltonian (J7J- On 
the other hand, to describe ID and 2D systems, the two- 
level Hamiltonian (J7J must be modified. The effects of 
reducing the number of spatial dimensions are twofold: 
the allowed values of the quantum number to become 
restricted, to to = for ID and to to = ±1 for 2D; and 
the interaction and tunneling energies are also modified. 
The 2D problem in fact has less symmetry than the 3D 
problem, since there is no axis along which the system is 
invariant under rotations. 



E. Characteristic parameters of the double well 
potential 

As a particular example of a ID double well potential, 
we will consider a Duffing potential 



V(z) = V Q 



8 



-z* + 



16 



2 4 + l 
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with barrier height Vo- The Duffing equation arises nat- 
urally in MF approaches to the double well problem, 
where chaotic oscillations of the atomic population at 
each well are found for time-dependent potentials [57ll58j |: 
the Duffing potential is also used in elementary portray- 
als of symmetry-breaking and phase transitions. We use 
this particular form of the potential to illustrate our re- 
sults. The results given in Sec. IIIII are valid for arbitrary 
double well potentials, since they depend only on the 
form of the Hamiltonian. Also, the criteria introduced in 
Sec. II VI hold for arbitrary double well potentials provided 
that the relevant hopping and interactions coefficients are 
properly calculated. Instead, in general, it is not possi- 
ble to represent these criteria for arbitrary potentials in 
a plane, as we do in Sec. [V] for the Duffing potential. 
Moreover, this potential gives a straightforward expres- 
sion for the double well potentials found in experiments 
and permits us to characterize the problem using, besides 
the number of atoms N, only two additional parameters. 



We show in the following that these two parameters arc 
related to the barrier height Vo, the distance between 
minima z m ; n , and the coupling constant g. 

For the Duffing potential © 32Vo/z^ in and, for 
lu x = u! y = u>, we have V x (x) = 16Vo x /z^ and likewise 
for V y . For an atom species of mass M, we numerically 
calculate the eigenfunctions tp n i m (r — rj), the energy lev- 
els Ei, and the hopping coefficients (|7d|) . Finally, for a 
number of atoms N, we completely characterize the sys- 
tem once g is known and the coefficients (|7ep are found. 
As we will see, we can distinguish different regimes in 
terms of criteria that relate the level spacing, the number 
of atoms, and the hopping and interaction coefficients. In 
this section, we describe a procedure to express these co- 
efficients in such a manner as to clearly identify these 
regimes for any atomic species. 

The recoil energy associated with a ID periodic 
optical lattice of wavelength A is defined as E r = 
2H 2 n 2 /MX 2 . Analogously, for the ID Duffing poten- 
tial, we consider A = 2z m j n and E r = h 2 n 2 /2z^ nin M. 
Dividing the potential by E r we can write V(z) = 
Vo {Sz 2 /zl in + 16z 4 /4i„ + 1) with 



Vo 
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Similarly, V x (x) = 16Vox 2 / z^ in , and likewise for V y (y). 
Then, 
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E r 
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where V(f) = V x (x) + %{y) + V z (z), (x,y,z) 
{l/z miTl ){x,y,z), and $# m (r) = z^ n ip im (r ± a) . Also, 

-g / dr|^ ro (r)| a |*^ m /(f)| a , 



E r 



with 



2E r z 
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min 



(10) 



Finally, the functions <frj£ m (r) and their corresponding 
eigenvalues can be numerically calculated as detailed in 
App. O Therefore, given the number of atoms N, we 
can characterize the problem with two parameters, i.e., 
Vo and g. Notice that the distance between wells z m i n is 
enclosed in the scaling procedure. 



F. Energy scales in the high and low barrier limit 

Let us now obtain, in the high barrier limit, the rela- 
tionships between the different relevant energy scales of 
the problem, i.e., the level spacing AE, Up!^,, and Ji m . 
In App. [B] it is shown, approximating the eigenfunctions 
with the spherical harmonics, that the interaction energy 
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of two atoms occupying the lowest energy level of one well 
is 



U = U™ = g 



(£)'■ 



[11) 



where a = Ak\/Vq. Also, the following relationships can 
be obtained among the hopping and interaction energies, 
in the high barrier limit: 

XJll = = Ultl = 2U& = (2/3)[/ 1 1 ° = U /2. (12) 

Finally, it is also shown in App. lAlthat the tunneling 
energies satisfy 



Jo 



(3/4)J] 



i±i, 



Jio 3> Jo 



(13) 



provided that ATS is not small. For notational simplicity 
we define 



Jo = Joo, and Ji = Ji, 



(14) 



Thus, in this high-barrier regime, there are four rele- 
vant energies: the energy level spacing, AE; the tunnel- 
ing energies, Jo and J\; and the interaction energy, Uq. 
Nevertheless, all these energy scales are not independent 
parameters since they depend solely on Vq and g. 

These results hold as far as the spherical harmonics are 
a good approximation of the eigenfunctions of the double 
well, i.e., for high Vq. For the low barrier limit, we use the 
expressions given in App. |Bj with numerical evaluation 
of the corresponding eigenfunctions and eigenvalues. In 
both limits, once the parameters, Vq, and g are given, all 
coefficients can be calculated, and, for a given number 
of atoms TV, all many-body eigenstates can be found. 
Thus the Duffing potential makes this a three parameter 
problem. Although we restrict our discussion to repulsive 
interactions, Uq > 0, our results hold for Uq < as well. 



G. Fock basis and dimension of Hilbert space 

Throughout our discussion, we operate in Fock space. 
An arbitrary state vector |$) in Fock space has the fol- 
lowing representation, 



Q-l 

l*> = Ci|*>F, |»)f 



i=0 



where 



n W i 

n jtm- 



;t \ 1 



,<*> 



|0>. 



(15) 



(16) 



Here VI is the dimension of the Hilbert space {|«)f}, i 
is the Fock-space index, and \ci\ 2 is the probability of 
finding n-J m atoms in the fth energy level of the jth 
well with z-component of angular momentum m when 
the system is described by state \ty). We work in the 



canonical ensemble, i.e., we require the total number of 
atoms 



(17) 



to be constant. Under this restriction, the dimension of 
the Hilbert space is given by 



fi = 



(N 



-1)! 



N\ (n m -iy. 



(18) 



where n m is the number of modes used to expand the field 
operator. For the double well in 3D we have £ truncated 
at 1, so n m = 8. For a large number of atoms, f2 scales 
like N 7 . 

The index i is chosen to increase with the number of 
atoms in well j = 1 of the lower level, with the number of 
atoms in the same well in the excited level with m = — 1 , 
then with m = 0, and finally with m = 1. Therefore, 
for the first N + 1 Fock vectors i = 1 + moo and they 
correspond to vectors with no occupation of the excited 
level. Thus, they satisfy 
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A + \ "100 /_ , \ "200 , , , 

(b\oo) (4oo) |0>, (19) 



for i = 0, 1, . . . , N . The one- level approximation 
can easily be recovered from Equation ([7|) by requiring 
i < N + 1. In this truncated space, the dimension of the 
Hilbert space reduces to that of the one-level approx- 
imation, namely, N + 1, the two- level Hamiltonian H 
reduces to the one-level Hamiltonian Hq, and we recover 
the LMG Hamiltonian. 

For the next 6iV vectors, one atom occupies the excited 
level and 



i =2 + nioo + N 



E 



(27V; 



10 



47V- 



li 



1) 



where Ne m is the number of atoms at level £ with z- 
componcnt of the angular momentum m. These vectors 
correspond to all combinations of N—l atoms in the lower 
level and a single atom in the excited level with m = 
— 1, 0, 1, in two wells. The Fock index i increases further 
with all combinations of p = 2, . . . , TV atoms occupying 
the excited levels and N — p atoms in the lower level. 



III. CHARACTERIZATION OF EIGENSTATES 

We begin our analysis with a characterization of the 
eigenstates \<f>k) of the two-level Hamiltonian ([7]). The 
eigenstates satisfy 



H\<f> k )=e k \tM, 



(20) 



where e k is the energy eigenvalue corresponding to the 
state \4>k)- The eigenstate label k is chosen to increase 
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with £fe. In order to describe these states, we will use the 
Fock-space amplitudes 



=i fc) = F <i|&>. 



(21) 



Insofar as the interlevel effects are not relevant, the 
eigenstates fall into one of two categories: harmonic- 
oscillator-like states (HO states) or MS states. When 
the barrier between wells is low, (/N > 1, all states are 
harmonic oscillator-like. This regime is known as the 
Josephson regime. On the other hand, MS states domi- 
nate the spectrum in the high barrier limit, f <C 1. This 
limit is known as the Fock regime. We recall that our 
coefficients U[™ and hopping coefficients, Jg depend on 
£ and to. Indeed, the tunneling coefficient for the ex- 
cited atoms with to = 0, J\, is much bigger than Jq. 
Then, between these two regimes, an mixed one can be 
found, for which MS states occur for atoms in the bottom 
level while HO states occur for the excited ones. When 
the level spacing is comparable to N 2 U~o, interlevel ef- 
fects can be no longer neglected, and another category 
of eigenstates emerges. These states show coupling be- 
tween MS states with atoms only in the lowest energy 
level and states with atoms in the excited one; we name 
them shadows of the MS states. 



and 



E 



Skin 



(24) 



where = n$ m + n$ m . Here, $ m is the one- 

level amplitudes and energies etmk are the one-level en- 
ergies. Both quantities can be obtained exactly in the 
non-interacting limit. The amplitudes are given by 



a K {k)[ 



(^) ^(n^N^pin^N^) , (25) 



where p(n^ m \N^) is the square root of the binomial 



'Kml 1 * Im i 

distribution, h„(k) (n^j^lN* 



im ' 



is a K^th order discrete 

Hermite polynomial, and o> K w(N^) is a normalization 

factor (for an expression of these coefficients see App. [C]) . 
Notice that the one-level Hamiltonians, expressed in the 
Fock basis, resemble, in the non-interacting limit, a har- 
monic oscillator potential truncated at hard walls. This 
gives rise to the binomial distribution and Hermite poly- 
nomial, as appropriate for such a potential. The corre- 
sponding eigenvalues are 



£klr, 



-Jim ( N. 



r(fe) 



Im 



(26) 



Non-interacting limit: harmonic-oscillator like 

states 



We first consider the limit C,/N — > oo by setting NUq = 
0, i.e., the non-interacting limit. In this case, both the 
energy levels and the orbital states are completely de- 
coupled. The two-level Hamiltonian Eq. is thus re- 
ducible to four independent one-level Hamiltonians, since 
all couplings between different £ and to depend on inter- 
actions. Furthermore, because [if, nu m + hiim\ = for 
all i and to, the eigenstates of the two-level Hamilto- 
nian must have definite occupation of the mth orbital 

( k) 

state of the £th energy level. Let be the number 

of atoms occupying the mth orbital state of the £th en- 



(fe) 



N. 



ergy level for the kth. eigenstate. Then, J^e m 
Let us denote the one-level eigenstates as \(j) K w), where 



• ' ^im * s tne one -l eve l eigenstate label. 
Then, the fcth eigenstate \cf>k) is a direct product of these 
one-level eigenstates: 



\4>k 



K 



(*0/ 



(22) 



Likewise, the ith Fock space amplitude c\ and the 
eigenenergy Ek can be expressed in terms of the one-level 
amplitudes and energies as 



i*0 



n 

£,m 



C iim <> N ( 



N 



(23) 



(k) 

Because the amplitudes c ilm resemble the eigenfunctions 
of the ID harmonic oscillator potential and the eigenval- 
ues £kim are linear in K$ , the eigenstates (|22[) are said 
to be harmonic-oscillator-like. 

The ground state \<po) is a coherent superposition of 
atoms in the lowest energy level of the left and right 
wells. The probability density of the ground state is 



,(o)i 2 



1 



,W I ( AT _ „W \ \ 

noo- 'noo J ■ 



with corresponding energy 

£ = -NJ . 



(27) 



(28) 



This result is readily generalizablc to an arbitrary num- 
ber of energy levels. 



B. High barrier: macroscopic superposition states 

We now turn our attention to the opposite, high 
barrier limit or Fock regime, C ^ 1- Let us as- 
sume first that Jo = Ji = 0, i.e., the infinite-barrier 
limit. In this regime Eq. (fT2"|) holds, and it is evident 
that none of the coefficients t//,™, can be neglected. 
Then, the eigenvectors of Hamiltonian ([7]) are not Fock 



vectors, due to the terms 



and 



J jlm 



(&jio) ^ii^-i-i +h- 
Nevertheless, we can neglect 
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these terms whenever 2AE ^> N 2 Uq. We will justify this 
criterion in the next section. Then, the eigenvectors of 
the resulting Hamiltonian are, indeed, Fock states with 
eigenvalues: 



(fc) 



/ in 



U 

m 1 



N, 



(fc) [ N e 



(A:) 



j,m" 



(29) 



up to (iV — 2i/)th order in Jq/Uq- In agreement with 
the rotational symmetries of the potential, the states 
l^i+i!^) and \4>^i,v) are degenerate. On the other 
hand, the energy difference between states 



4 (±). 
u oo > 



v) and 

l^^; v) is on the order of NAE when NU° < 2AE. 

Since J\ > Jo, it is also possible that Jo <§C Uq but 
Ji > Uq. Then, atoms in the bottom level behave as in 
the Josephson regime while the ones in the excited level 
behave as in the Fock regime. We will show numerical 
examples of these mixed regime in Sec. [V] 

C. Shadows of macroscopic superposition states 



According to Eq. (p?9")) the number of degenerate eigen- 
states depends on the occupation of the excited level. 
For example, for no atoms in the excited level, there 
are two degenerate eigenstates obeying n[ k Q = n£2 and 



,(fc) 

'200 



'100 • 



'200 



For one atom in the excited level with 
m = there are four degenerate eigenstates, since the 
excited atom can be located in any of the two wells, thus 
giving four combinations. Let us now consider Jo <C Uq 
and J\ -C Uq. Non-degenerate perturbation theory gives, 
in every case, that the eigenvectors are quasi-degenerate 
symmetric and antisymmetric combinations of the corre- 
sponding Fock vectors (see App. [D]). For the particular 
cases in which all atoms occupy the same level, and for 
which the angular momentum of each atom is oriented 
along the z-axis, that is, m = ±£, the eigenstates are 



\0 



(30) 



for < v < N/2; we have neglected terms on the order 
of {Jn/Ua)*- 2 " and smaller. Here 



V2 



1 



u Um ) [ u 2hn 



± 



y/v\(N -V)\ 



N—u 



|0>, 



(31) 



Let us turn now our attention to the effects of the cou- 
pling between energy levels in the high barrier limit or 
Fock regime, Jo -C Uq and Ji -C Uq. Let us consider the 
hopping terms and the interaction terms that account 
for same-site interlevel hopping in Hamiltonian ([7]) as a 
perturbation to the decoupled Hamiltonian. As shown 
in App. [E] for the N + 1 eigenvectors with zero occupa- 
tion of the excited level, the first order approximation 
to the eigenvector shows coupling to states with differ- 
ent number of atoms in the excited band. These cou- 
plings are associated with the destruction of two atoms 
in the lower level and creation of two atoms in the first 
level with m = or one with m = 1 and the other 
with to = — 1. The corresponding coefficients are neg- 
ligible as far as 2AE > NU . Nevertheless, if 2AE is 
comparable to NUq , Fock vectors with non-zero occupa- 
tion of the excited levels are coupled to the MS states. 
We call these coupled vectors shadows of the MS states 
\4>± ; Hi). Similar results hold for MS states with nonzero 
occupation of the excited level. Coupling between differ- 
ent levels in asymmetric double wells or optical lattices 
plays a fundamental role in far-from-equilibrium dynam- 
ics showing Landau-Zener (LZ) coupling [5^,[6^|. Here 
we focus on statics and on the symmetric case, leaving for 
future work the study of how LZ coupling between differ- 
ent wells in asymmetric double well potentials is modified 
in this regime. 



is an MS state in which v and N — v atoms simulta- 
neously occupy the mth orbital state of the £th energy 
level of both wells. Here ipo is the usual arbitrary phase 
associated with vectors in a Hilbert space. We will set 
ipo = for the rest of this Article. The special case v = 
represents an extreme MS state in which all N atoms si- 
multaneously occupy the left and right wells. These MS 
states can be either symmetric (+) or antisymmetric (— ). 

The eigenstates (f3"U)) occur in nearly degenerate pairs of 
symmetric and antisymmetric MS states. The level split- 
ting between the states |<^^;i>) and is Aet{v) 
where 



WttlJlmW'Z^-fjN -V)\ 

i/![(7V-2i/-l)!] 2 



(32) 



IV. BOUNDS ON THE USE OF A ONE- AND 
TWO-LEVEL APPROXIMATION 

We have distinguished two main regimes, the Joseph- 
son regime, in which the eigenstates arc HO-like and, 
the Fock regime, in which MS states can be found. The 
Josephson regime is characterized by: 



NU, 



tm 



< l. 



The Fock regime is characterized by: 

J, 



'/■Fock, 



jjlra 
u lrn 



< 1. 



(33a) 



(33b) 
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Notice that, as stated above, these criteria should be 
evaluated for both levels. Then, it is possible that the 
Fock regime holds for atoms in the bottom level, while 
the Josephson regime holds for atoms in the excited level. 
Since J\ > Jq the contrary is not true. Hence, in general, 
we distinguish three regimes: the Josephson regime, the 
Fock regime, and the mixed regime. In the first two 
regimes, the corresponding criterion holds for atoms in 
both levels. In the latter, the Fock criterion holds for 
atoms in the bottom level, while the Josephson criterion 
is satisfied for atoms in the excited level. Let us show, 
for these regimes, the bounds on the one- and two- level 
approximations. With our choice of indexing states, the 
one-level approximation corresponds to truncating the 
size of the Hilbert space to N + 1. Then, the bounds 
we present below will describe the regime in which this 
truncation is valid. 

Let us consider first the Josephson regime. The energy 
levels are coupled by the interaction energy Up^, t £ ^ £' . 
In this regime all coefficients Up 1 ^, are small, and then, 
the coupling between levels is weak. Energy levels only 
become completely decoupled when £/|/™ , = 0. However, 
let us show that eigenvalue crossings arc induced by the 
presence of the excited level. Let us assume that Uq = 0, 
which implies that U^, = 0. According to Eq. (J25J), the 
maximum of the eigenvalues for the first N+l eigcnstates 
with no occupation of the excited level coincides with the 
minimum of the eigenvalues of the states with one atom 
in the excited level, if 



AE 



XJc 



Jo (2N - 1) + J l 



1. 



(34) 



this being the criterion that determines the first eigen- 
value crossing in this regime. For xj os > 1 no crossing 
occurs. Moreover, for 



os,gs 



AE 

Jl — Jq 



1, 



(35) 



the first crossing involving the ground state occurs, i.e., 
the ground state shows non-zero occupation of the ex- 
cited level if XJos, gs < 1- 

Analogously, in the Fock regime, the first eigenvalue 
crossing occurs when the condition 



XFock 



'00 



(A 2 + 2A - 3) 



2(AE- 



2U™(N-1)) 



= 1 



(36) 



is met. This condition is obtained equating the maximum 
eigenvalue given by Eq. (j2"9")l for states with no occupation 
of the excited level, to the minimum eigenvalue given by 
this equations for states with one atom in the excited 
level. For XFock < 1 no crossing takes place. For large N 
this criterion can be approximated by N 2 U~o ~ 2AE. On 
the other hand, the first crossing involving the ground 
state occurs when the condition: 



XFock, gs 



(3/2)U™(N- 



1 



2U™(N 



1) 



AE 



1 (37) 



is satisfied. For XFock.gs < 1 the ground state shows non- 
zero occupation of the excited level. For large N this 
criterion is NUq ~ 2AE. 

If we consider Eq. (|12p . in the high barrier limit, the 
criteria (|3^|) and (f3"T)) turns into: 



and 



approx _ Uq(N 2 — 1) _ 
X F ock - 2AE ~ 



approx 
XFock, r- 



2AE 



(38) 



(39) 



For completeness, let us consider the following crite- 
rion: 



Xshadow 



V2U&y/N(N-l) 



U (AN - 6) 



2f/ 10 



8U™(N 



2AE 
(40) 



In the Fock regime, the expression of the MS states in 
terms of the Fock basis gives only two relevant coeffi- 
cients, those corresponding to the two Fock vectors that 
are superimposed, as detailed in Sec. IIIII But, as the in- 
teractions are increased the coefficients corresponding to 
the shadows of MS states are more relevant. The great- 
est of these coefficients corresponds to a Fock vector with 
two atoms in the excited level with m = 0. The crite- 
rion (|4"U)) is the expression of this coefficient, which is 
obtained in App. [E] using perturbation theory. Hence, 
for small Xshadow we can neglect the interlevel coupling, 
since all coupling to other Fock vectors will be negligi- 
ble. If we consider the relations given by Eq. (TT2"j) . this 
criterion can be written as: 



approx 
Xshadow 



U y/N(N-1) 

V2AE (U /2AE - 2) ' 



(41) 



Moreover, the single particle eight-mode basis turns to 
not be appropriate to express the field operators if the 
barrier height is smaller than the energy gap between 
levels AE, i.e., if 



XAE 



Vq 
AE 



< 1. 



Finally, the criterion 



Xweak — 



N*U 
2AE 



< 1 



(42) 



(43) 



has to be met to account for the weakly interacting gas 
condition, Eq. ([2]). Equation (|43|) is obtained using the 
analytical form of Uq and AE to elliminate ah /a s in 
Eq. ([2J, taking into account all the scalings performed. 

Therefore, the ground state shows occupation of the 
excited level when criteria (ptS"j) and (|3"T|) . or its high bar- 
rier version (|39[) . are fulfilled. Moreover, we are interested 
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in describing cat-like MS states which are typically ex- 
cited eigenstates. Therefore, a characterization of eigen- 
value crossings of energies other than the ground state is 
also relevant. These crossings appear when criteria (|3~4")l 
and (|36|) . or ([38]) , are satisfied. Also, criterion (|40|) . or 
(|4l"jl . indicates the presence of shadows of MS states. 
It is important to notice that, for large N, both crite- 
ria d37J) and (gOl) becomes NU/2AE, while criterion ((55)) 
becomes N 2 U/2AE. Then, concerning MS states, the 
excited level plays a relevant role insofar N 2 U ~ 2AE. 
Furthermore, when NU ~ 2Ai? the ground state shows 
occupation of excited levels and the coupling between 
levels is non- negligible. Then, for large N, the use of 
single-particle wavefunctions and only a few energy lev- 
els is appropriate to the regime 



NU, 



X mo del — 



2AE 



< 1. 



(44) 



If this condition is not met, our approximation is inac- 
curate and alternative treatments become necessary 
l32j . Notice that this criterion is identical to the regime 
for which mean field theory is valid. Finally, criteria (|42|) 
and (|43[1 are also two limiting criteria for the model. 



V. NUMERICAL RESULTS 

Let us use the previous criteria to completely demar- 
cate the different regimes in the Vo-g plane. At every 
point of this plane different values of the interaction co- 
efficients, Ufj^,, the hopping coefficients, Jg m , as well 
as the energy level spacing, AE, are obtained. Since 
all criteria depend on these parameters, we can draw in 
this plane the curves for which the different criteria are 
satisfied, thus determining different regions in which the 
eigenstates have been characterized. We will illustrate 
the results with examples obtained after exact diagonal- 
ization of the Hamiltonian ([7]). Also, we will determine 
in that plane the limits of validity of our model. 

In figure[2]we represent, for N = 8 atoms, the curves in 
the Vg-g plane for which the criteria that characterize the 
first crossing, Eq. and Eq. (|3"6"|) . are fulfilled. These 
criteria are valid for large g and Vo, respectively. On 
the other hand, an intermediate regime arises between 
the Josephson and Fock regimes for which the interac- 
tion coefficient is comparable to the hopping coefficient. 
Since the results obtained theoretically are not valid in 
this intermediate regime we have performed exact diago- 
nalization of the Hamiltonian ([7]) to determine the pairs 
Vo-g that lead to the first crossing. The corresponding 
interpolated curve is presented as the criterion Xnumeric 
in the figure. In Fig. G2 we also represent this curve. 
Points 1, 2, and 3 in this figure correspond to Figs. 2] (a), 

(b), and (c), where the coefficients |c^J 2 for the first 
N + 1 + 6N eigenvectors are represented, thus showing 
the first crossing in the first N + 1 eigenvectors. Hence, 
there is no crossing in the region to the left of this curve. 
We also represent in Fig. [5] the curve for which criteria 




FIG. 2: (Color online) Characterization of the crossings and 
validity of the model in the Vo-g plane (both dimensionless) , 
for the Duffing potential. Different criteria that characterize 
the first crossing and the first crossing involving the ground 
state, as well as the limits of validity of the model for N = 8 
atoms. Three main areas are distinguished: (i) the area in 
which no crossing occurs; (ii) the area in which eigenstates 
with occupation of the excited level emerge among the first 
jV + 1 eigenvectors; (iii) and the area in which even the ground 
state shows occupation of the excited level. The one-level 




-2 

iog 10 (g) 

FIG. 3: (Color online) Characterization of the different 
regimes in the Vo-g plane (both dimensionless) . Different cri- 
teria that characterize the Josephson and Fock regimes for 
both levels when N = 8. The numbered points represent the 
examples given in subsequent figures. For completeness, the 
three regions distinguished in Fig. [2] are also represented. 



for the first crossing involving the ground state in the 
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FIG. 4: (Color online) Eigenstates for three examples of the 
first crossing in the Josephson, intermediate, and Fock regime. 
Probability amplitudes |c'^J 2 of the first N + 1 + 6N eigen- 
states for (a) Vo = 20; g « 1CT 4 , (b) V = 20.6; g^2x 10~ 3 , 
and (c) V = 49; g « 10~ 2 , when JV = 8 atoms. Above 
the dashed line one atom occupies the excited level. Dot- 
ted lines separate the Fock vectors for which this atom shows 
m = —1, 0, 1, respectively. 



Josephson [Eq. (|35j) ] and Fock [Eq. (l37l) ] regimes are sat- 
isfied. In the region between these curves and the one 
given by Xmimeric the crossings do not involve the ground 
state. Finally, the curves associated with criterion (|4"3")l 
and criterion (|42]) demarcate the limits of validity of the 
model. 

Once these three main regions have been defined, let us 
show in this plane the different regimes in which the indi- 
vidual eigenstates have been characterized. We represent 
in Fig. |3j for the bottom level and for the excited level 
with m = 0, the criterion for the limit of the Josephson 
regime, Eq. (|33a[) . for t7j s. ^ = 0.1, i = 0, 1. Similarly, we 
represent the criterion for the limit of the Fock regime, 
Eq. ()33b[) . for f?Fock,£ = 0.1, £ = 0,1. Cusps in both 
curves are an artifact associated with the resolution of 
the interpolation between the points in which we have 
numerically calculated the eigenfunctions and all the co- 
efficients. 

To further characterize these regimes, we use the eigen- 
values of the single particle density matrix, whose ele- 
ments are defined as: 

Pij = (*o|S&|*o>, 

where <Fo is the ground state. Notice that, since the 
indices j£m permit one to run over all combinations j, £, 
and m, this matrix is of dimension eight. The curve Ai 
represents the pairs Vo-g for which the largest eigenvalue 
takes the value Ai = 0.85. Also, the second eigenvalue 
starts to grow to the right of this curve. Thus, this curve 
defines the limit of the Josephson regime for the bottom 
level. In Fig. [5] (a), (b), and (c) we represent examples of 
the coefficients \c^L | 2 for the first N+1+6N Fock vectors 
and eigenstates in the Josephson, intermediate, and Fock 
regimes, respectively. These examples are represented in 
Fig. [3] as points 4, 5, and 6, respectively. 

In Fig. [6] we show that the example corresponding to 
point 7 in Fig. [3] belongs to the mixed regime. Hence, the 
first N+l eigenstates, which show occupation of only the 
bottom level, are MS states. Conversely, all eigenstates 
corresponding to occupation of only the excited level with 
m = behave as HO states. Notice that the latter are not 
consecutive eigenstates, since there are many crossings in 
the excited level that we have not considered in detail. 

On the other hand, in the Fock regime, as g is in- 
creased shadows of the cat states may appear, as dis- 
cussed above. The curve for which criterion f4T)]) is sat- 
isfied for Xshadow = 0.15 is represented in Fig. [3J In 



Fig. [7] (a) and (c) the probability amplitudes |i 



(fc) 12 

iim I 

for two examples showing shadows of MS states are 
shown. These two examples are represented in Fig. [3] 
as points 8 and 9, respectively. Fig. [7[b), represents 
the coefficients for the eigenstate with index k = 120 
of the first example. This state is a superposition of 
the Fock vector |0, 8) £3;=i, TO |0,0), with Fock index 
i = 1 and the Fock vectors |6, 6) <g> |0, 1) <g> |0, 0) <g> |0, 1) , 
|0,6) <g> |0,0) ® |0,2) <g> |0,0) , with indices i = 107 and 
i = 135, respectively. These vectors are also coupled to 
the Fock vectors with indices i — 9, i — 134 and i = 155, 
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FIG. 5: (Color online) Eigenstates for three typical exam- 
ples of the eigenvectors from the Josephson, intermediate, 
and Fock regime. Probability amplitudes c^J 2 for the first 
TV + 1 + 6TV eigenstates for TV = 8 atoms, Vo = 32.4, and (a) 
g w 1CT 4 , (b) g « 8 x 1(T 4 , and (c) g ^ 7 x 1CT 3 . Above 
the dashed line one atom occupies the excited level. Dot- 
ted lines separate the Fock vectors for which this atom shows 
m = — 1, 0, 1, respectively. 
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FIG. 6: (Color online) Eigenstates for the mixed regime. 
Probability amplitudes |c^J 2 (a) for the first TV + 1 eigen- 
states, in which only the bottom level is occupied and (b) for 
all eigenstates corresponding to occupation only of the excited 
level with m = when TV = 8 atoms, Vo = 37, and g ss 10 -3 . 
Notice that, while the former eigenstates are MS states, the 
latter are HO-like states. 



corresponding to the MS of the previous ones and those 
with the same number of atoms in the other well. As g 
is increased further, the shadows of MS states are more 
relevant, being noticeable also for the ground state. In 
Fig. [Hd) we represent the coefficients for the ground 
state of the second example. In this case the ground 
state displays non-negligible coupling to its shadow Fock 
states. The ground state is then the superposition of the 
Fock vector |4, 4) <E)i=i. m [0, 0) , with Fock index i = 5, 
to Fock vectors with indices i = 111, i = 130, i = 139 
and i = 151, i.e, the vectors |4, 2) ® |0, 1) ® |0, 0) <gi |0, 1) , 
|4, 2) <g> |0, 0) <g> |0, 2)® |0, 0) , and the similar ones obtained 
after interchanging the well index. Notice that if g is in- 
creased further crossings involving the ground state will 
take place. 

Finally, let us show how this scenario changes as TV is 
increased. In the Fock regime, the curves that represent 
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FIG. 7: (Color online) Shadows of excited MS states and of the ground state, (a) Probability amplitudes |c'^J 2 for the first 
(N + 1) + 6iV + 21(N — 1) Fock vectors and eigenstates when Vo = 38.4, j~3x 10 -2 , and N — 8 atoms. One atom occupies 
the excited level above the dashed line at i = N + 1. Two atoms occupy the excited level above the second dashed line at 
i = N +1 + 6N . (b) Coefficients for the 120-th excited state, the most excited state with all atoms in the bottom level. This MS 
state shows non- negligible coupling to Fock vectors with two atoms in the excited level, which are its shadows, (c) Probability 
amplitudes |c^J 2 for the first (N + 1) + 6N + 21(N - 1) Fock vectors and eigenstates when Vo = 49.8, g « 7 x 10" 2 , and 
N = 8. (d) Coefficients for the ground state for this case. Here, even the ground states shows non-negligible coupling to Fock 
vectors with two atoms in the excited level, being these its shadows. 
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criteria (|57| . ([ID]), and (|43p. each tend to a straight 
line as Vq is increased. We have checked numerically that 
as N is increased these curves move to the left without 
changing their slope. In Fig. EJa) we represent how the 
value of g for Vo = 100 changes as N is increased. Then, 
the area for which the one level approximation is valid 
is reduced as N grows, while the region for which shad- 
ows of MS states can be observed is increased. Using 
the harmonic oscillator approximation, it is possible to 
obtain expressions for the criteria (J3BJ, (|3T|) . and (|4"3")l : 

approx 1 approx 

*F°ck ~ p(N 2 - 1)' X Fock, s - p (JV- 1)' 



and 

approx 

X weak ~ pN 1 / 3 ' 

where p — tt 3 ^ 2 Vq^ 4 /\/32. For criterion an expres- 
sion cannot be obtained, but a recurrence relation is ob- 
tained. These expressions have been validated numeri- 
cally, showing good agreement with the numerical curves. 

On the other hand, we have checked numerically that 
the curve that represent criterion (f34|) . valid for the 
Josephson regime, moves upward without changing its 
slope as N is increased (notice that the other two do not 
depend on N). In Fig. [SJb) we represent how the value 
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of Vo for g = 10 -7 changes as N is increased. Once more, 
the area for which the one level approximation is valid is 
reduced as N grows. 



VI. SUMMARY AND DISCUSSION 

We developed a Fock space picture of the stationary 
states of a system of ultracold bosons in a 3D double- 
well potential using a two-level, eight-mode approxima- 
tion. These modes are 3D single particle eigenfunctions 
with on-well angular momentum I and z-component of 
the angular momentum m. We have identified all the 
processes relevant in such a picture. First, familiar pro- 
cesses occur, such as interaction of atoms in the same 
well and the same or different level, or hopping betwen 
atoms in the same level and different wells. On the other 
hand, other less common processes play a fundamental 
role. These are the hopping between pairs of atoms in 
the bottom level with m = and atoms in the excited 
level with m = and hopping between pairs of atoms in 
the bottom or excited level with m = and one atom in 
the excited level with m = 1 and other with m = — 1. We 
showed that these hopping processes are related to the 
interaction and not the hopping coefficients. Therefore, 
in addition to the level spacing, AE, and the hopping 
and interaction coefficients in the bottom levels, Jo and 
Uq, other coefficients have to be considered. These are 
the interaction coefficients between atoms with the same 
or different values of £ and m, U[™ , and hopping coef- 
ficient between atoms in different wells, Je m - 

We found that all the coefficients are closely related 
and, indeed, in the high barrier limit, they can be deter- 
mined in terms of only three of them, Uq = Uqq, Jo = Jqq, 
and J i = Jio- sec Eqs. (fT12j) and (fTUj) . Nevertheless, in 
the general case, the coefficients have to be evaluated 
numerically. They depend on the particular form of the 
double well potential and on the coupling constant g. In 
this Article, we chose a Duffing potential to illustrate our 
results numerically, which allows one to reduce this de- 
pendence of the coefficients on the particular geometrical 
form of the potential and the coupling constant to only 
two parameters, Vq and g, defined in Eqs. (0 and ([TU| . 
respectively. V) is related to the barrier height Vo and 
the distance between wells z m i n , while g is related to the 
coupling constant g and z m - m . 

BECs in double wells were previously thought to ex- 
hibit only two main physical regimes, the Josephson and 
the Fock regime. These are characterized in terms of the 
hopping coefficient, the interaction coefficient and the to- 
tal number of atoms N. Wc derived new relevant coef- 
ficients in the problem, and consequently a number of 
new regimes were identified. Although many new coef- 
ficients are considered, they are determined in terms of 
only two parameters, allowing one to consider the double 
well problem in terms of three parameters: Vq, g, and 
the number of atoms, N. For certain values of the coeffi- 
cients the eigenstates show a HO-likc behavior, while for 



others they are MS states corresponding to the conven- 
tional Josephson and Fock regimes, respectively. These 
regimes have to be considered separately for the bottom 
and the excited level. Moreover, we showed that MS 
states with non-zero occupation of the excited level can 
occur. These excited MS states show angular momentum 
degrees of freedom. Finally, we found a mixed regime, in 
which the eigenvectors with no occupation of the excited 
levels are MS states, while the eigenvectors showing solely 
occupation of the excited levels are HO-like states. We 
found criteria to distinguish all these regimes, which, for 
fixed N, were represented as areas in the Vo-g plane. 

The eigenvectors are also different in another region, in 
which coupling effects between levels become important 
since the interaction energy is comparable to the energy 
level spacing. In this regime, the interaction energy is 
much greater than the tunneling energy and the eigen- 
states are MS states which mix energy levels, showing 
shadows of cat-like states. Lowest order perturbation 
theory couples states with atoms solely in the bottom 
level to states with atoms in the excited level. We found 
the criterion to distinguish this region and represented 
the corresponding curve in the Vo-g plane. 

Moreover, eigenstates involving occupation of the ex- 
cited level can emerge among the lowest lying eigenstates, 
for certain values of the interaction and hopping coeffi- 
cients. We found the criterion that permits one to iden- 
tify whether such a crossing takes place for a given set 
of interaction and hopping coefficients. This criterion 
permits one to define the region in the plane Vo-g for 
which such a crossing does not occur. For certain values 
of the relevant coefficients, even the ground state of the 
problem shows occupation of the excited level. The cri- 
terion for such a ground state to exist was found. Again, 
this determines another region in the plane Vo-g in which 
this ground state with occupation of the excited level 
does not occur. Consequently, three main regions were 
identified: the region for which no state with occupa- 
tion of the excited level emerges among the first N + 1 
eigenstates, the one for which this occurs, and the re- 
gion for which even the ground state shows occupation 
of the excited level. For large N, the criterion that de- 
termines the first crossing for the excited MS states is 
approximated as N 2 Uo/2AE while the criterion that de- 
termines that the ground states shows occupation of the 
excited level is approximated as NUq/2AE, Also, the 
criterion obtained for the shadows of cat states to be rel- 
evant also can be approximated as NUo/2AE. Then for 
N 2 Uo ~ 2AE a two-level approach is necessary to study 
excited MS states while for NUq ~ 2AE our approxi- 
mation is inaccurate and alternative treatments become 
necessary [3ll. l32l| . 

Finally, we have established how all these criteria 
change with N. As the number of atoms is increased, 
the area for which the one-level approximation is valid is 
reduced. Also the region corresponding to the Josephson 
regime is reduced, showing, as expected, that mean field 
approaches are valid for small interaction Uq, provided 
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that N is not big. It may appear counterintuitive that 
MF approaches to the problem are less valid as N is in- 
creased, once Uq is fixed, but one must take into account 
that the criterion for the limit of the coherent or Joseph- 
son regime is NU/J -C 1, thus involving both variables. 
For completeness, the limits of validity of the model have 
been clearly identified and represented in the same Vg-g 
plane. 

The Fock picture for ultracold bosons in the 3D double 
well potentials developed in this Article will be used in 
the future to gain insight in the study of dynamical tun- 
neling ultracold bosons in 3D double wells. Now that all 
regimes have been clearly identified, one can study the 
tunneling properties of different initial population imbal- 
ances in the different regimes. Then, the possible initial 
states can show occupation of the excited level and it is 
expected that complicated and rich dynamics emerge out 
of the different regimes. 



where a = Att^/Vo. The corresponding eigenvalues are 
eo = 2a/n 2 and e\ = Aa/n 2 . The same eigenfunctions 
and eigenvalues are valid for the other two coordinates 
y and z. Then, the ground state of the single particle 
Hamiltonian (|A1|) is 



-ar 2 /2 



with eigenvalue Eq = 3eo = ■ Analogously, the first 
excited eigenfunctions are 



®io{r,0,<f) = <t)a{x)4>o{y)(j3i{z) 

3 

4 a?re- ar2/2 cos(9), 



\ 7T 



and 
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Appendix A: Eigenfunction of the single particle 
Hamiltonian 



The eigenfunctions $(f) of the single particle Hamil- 
tonian 



1 



V(r), 



(Al) 



permit one to obtain the coefficients Ji m /E r and 
UfT^,/E r after the scaling described in Sec. Ill El 
Since the potential is separable we can write $(f) = 
(fi(x)4>(y)<fi(z), where 4>{x), <j>(y), and (p(z), are the eigen- 
functions of corresponding one dimensional potential 
with eigenvalues e x , and e z respectively. 

As discussed in Sec. Ill CI for the high barrier limit 
$nim(r, 0, <p) = Rne(r)Ye m (6,(p). For convenience, let us 
write these functions in terms of the ID functions 
<j)(x), 4>{y), and <fi{z). The first two lowest excited ID 
eigenfunctions are 



and 



4>i (x) 



Ct 4 _ „=2 



x e 



ax 1 jl 



(A2) 



(A3) 



$i±i(r,0,<p) = -j=[<j)i(x)(i3o(y)(j)o{z) 



± i(j)o(x)(j)i(y)(j) (z) ^ 



with eigenvalue E\ = 2eq + ei = |§. Notice that we 
have included in the last expression the Condon-Shotlcy 
phase convention, as common in the quantum mechanical 
literature. 

For the low barrier limit this approximation is no 
longer valid. Nevertheless, we can proceed in the same 
manner to find the eigenfunctions and eigenvalues numer- 
ically. Let us consider the numerical solution localized at 

well j, calculated as $ (z) = ^= \ i>\{z) ± ^f(^)) > where 

ipi(z) and ^2(2) are tne nrs t two numerically calculated 
eigenfunction of the Duffing potential with eigevalues e" 
and £§, respectively. These functions have been calcu- 
lated with a imaginary time relaxation method. If we 
consider the plus sign j = 1 , while j = 2 in the other 

case. Similarly <j){(z) = ^3 (z) ± V>f(.2)) 7 where 

ip^(z) and ^}%{z) are the third and fourth numerically 
calculated eigenfunctions of the Duffing potential, with 
eigenvalues eg and ef, respectively. We can obtain $>ji± m 
using this numerical functions for the ID functions in 
the z variable in the expressions given above for the 
high barrier limit. The corresponding eigenvalues are 
£00 = 4a/TT 2 + (e? + eg)/2 , E w = Aa/n 2 + (eg + e%)/2, 
and^iii = 6a/7r 2 +(e?+ef)/2. Figure H shows an exam- 
ple of the numerically calculated localized eigenfunctions 
for the low barrier limit using this procedure and the ana- 
lytical approximation valid for the high barrier limit. As 
shown, the numerical eigenfunctions are deformed in the 
z direction compared to the analytical ones, thus giving 
higher values of the hopping coefficients Je m - 
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Appendix B: Expressions for the coefficients 

Let us use the expressions given in App.lAlto calculate 
the interaction coefficients. Accordingly, we have 

lC = g J d 3 f|$ 00 (f)| 4 = 5 C/ XC/ 5 , 

where U§ = J dx |</>o(i)| 4 , and similiarly for Uq and E/q. 
In the low barrier limit, we numerically evaluate 4^{z). In 
the high barrier, we use also the analytical approximation 
for all the ID functions and then 



TI 00 - n 
u oo — 9 



a \ 2 
2tt, 



Analogously, we find 

UlU = g I d 3 f|l> 1±1 (f)| 2 |$ 1±1 (f)| 2 



where TJ§\ = / rfi|</ , o(5)| 2 |0i(ai)| 2 | an( l Uf is definded as 
above. Similar expressions hold for Uq X and U^. Finally, 
wc have 



U$ = g I ^r|<M?r = WW 



U$ = g I d 3 i\^oo(i)\ 2 \^io(i)\ 2 = gUSU^, 



U, 



i±i 
oo 



g J rf 3 f|$oo(f)| 2 |$i±i(f)| 2 

{{ukulul + ululul 



uir = 9 d 3 fi$ 10 (f)i 2 i<i> 1±1 (f) 



— ± \TT X TT y TT z 4- TT X TT V TT Z 

— 2 ' 01 01 01 01 



In the high barrier limit, the previous expressions give 
Eq. (|12j) . In the low barrier limit, we numerically evaluate 
4>i(z). 

On the other hand, the hopping coefficients are 

-^v 2 + y(f) ) $, 



E Q I d 3 f$| 00 ^v 00 ) = E J§, 



where Jq = J dztpQ* (jy^ . In the high barrier limit, 
6 e~ Q / 4 (l2 + a(a - 4)V - ( °~ 2) ° 3 ) 



'00 



Similarly 

Jio = 



— V z + V(r) $ 



where Jf = J dz(ft[*(fr{ . In the high barrier limit, 

4 e -«/4 



J 



10 



■{a[36 + a(a-6)]Vb 



[l2-a(a-12)]a 3 



120} 



Finally, 



J 1±1 = J dr$* 1±1 ( --y 2 + V(r) ) <F /1± i 



According to the expressions for the eigenvalues given in 
App.O in the high barrier limit Ji±i = 4/3 Jqo- 



Appendix C: Low Barrier Limit 



The coefficients c^) m defined in Eq. (|2"3"]l . are associ- 
ated with the zth eigenstate, with atoms in level £, 
^-component of the angular momentum to, nfl atoms 
in well j = 1, and — nfl in well j = 2. Equa- 
tion (|25p gives their expression in terms of the the bino- 
mial coefficient p(n^ m \N^) and the normalization con- 
stant a K w(NfrH), which arc 



V (r 



II n 



N 



1<2n\ n (i) UN {k) -n {i) V 



and 



(N^ _ /r (fe) V 
\ iV„ (fe) 'id fe) l 



respectively. 



Appendix D: High Barrier Limit 

The unperturbed Hamiltonian is 



j,m I. rn' 



(Dl) 
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The cigenfunctions of this Hamiltonian are the Fock basis 
vectors, with eigenvalues given by Eq. (|29p . The perturb- 
ing Hamiltonian is Hj = - J2e, m Hj,e m , with 

Hj,£m = — ^j>2 J (m E bjerrfij'&m- 
j,m jjtj' 

The dimension of the degenerate subspaces, while de- 
pending on the number of atoms in the excited level, 
is always a multiple of 2. The perturbing Hamiltonian 
Hjj m acts on these 2-dimcnsional subspaces. Hence, we 
can diagonalize it in each subspace. Since 



P 





Vk'\Hj,£ m \<Pk 



Pk\Hj,i m \(p k > 




with 



for s < N, 



(k) 



<Pk\H'jjJ<p*)=O t 



2n uL we must use Mm 



order degenerate theory. In such a case, we obtain 



P' 




with 



P' - 

M2 — 



and 



(k) 

nim^N^ - n e 



H 



j,im 



N, 



(k) 



ne in , n£ m 



(0) _ (fc)(0) 
t nie m 



P' 



2Ui m ( N {k) - n [k) \ ' ( Jtm \ 



Nr'-2n\ 



"•Urn- 



N (k) _ 2 (fc) 



Therefore, the eigenvectors and eigenvalues of the matrix 
are 



and ±P{ 2 . 



The eigenstates include superposition states of the form 




"l!in: iv ta "-Urn 



± 



N, 



(fc) 
tin 



„(*0 ( fe ) 



and consequently symmetric and antisymmetric states 
appear in nearly degenerate pairs with energy differences 
equal to P[ 2 - The eigenstate of the complete Hamiltonian 
is the direct product of these states, with eigenvalue equal 
given by Eq. (|29p . plus (minus) the energy differences P[ 2 
for the symmetric (antisymmetric) state considered. 



Appendix E: Interlevel perturbation theory 



We consider the unperturbed Hamiltonian as 



) j,m" v 

and the perturbing one is then 



■j£m' ^mm') (2 ^mm') 



n j00 Tljlm" > , 



6 ]io) hn*>ji-i + h - c - 



Hp = J2\ (! - S mo)UH 

j,m I 

+ 1 U ™' (*Joo) W' hl-m" + h 



j,m" 



Let us illustrate the interlevel coupling due to the per- 
turbing Hamiltonian for the N + 1 eigenvectors with zero 
occupation of the excited level, i.e., states of the form 
|0(°>) = \^;n W o)<S, jihm \^;0). The first order ap- 
proximation to the ground state \4>^) gives 



|0< 



where 



C 2 



C3 



C4 



: ci |moo - 2, 7i 2 oo) O |0, 0) ® |2, 0) g> |0, 0) 
c 2 |nioo, moo - 2) g> |0, 0) ® |0, 2) ® |0, 0) 
c 3 |moo - 2, n 200 ) g> |1, 0) ® |0, 0) ® |1, 0) 
c 4 |nioo, moo - 2) g> |0, 1) ® |0, 0) ® |0, 1) , 



Coo a/2 (moo) (moo -1) 



LT™ (6 - 4n 



100 J 



2^ 



8 J/™ + 2AE' 



U$ ^2(71200)^200 - 1) 



(6 - 47i 200 ) + 2C7 1 1 ° + 8E7™ 

rl±l 



2AE' 



Cqq y^lOo) ("-100 - 1) 



4nioo) + 8L/ 1 ° 1 + 2A^' 



C Qo 1 ^(^200) (™200 - 1) 



u; 



00 ( 6 - 4 "2 



2AE 



Notice that coefficients C3 and C4 are larger than the other 
two. Then, the most relevant modification in the eigen- 
vectors is the coupling to the vectors jnioo — 2, 71200) ® 
|l,0)<g)|0,0)(g)|l,0) and |ni 00 , n 2 oo - 2) ® |0, 1) <8 |0, 0) ® 
1 0, 1) . For large N the expressions of these coefficients 
become NUq/2AE. Therefore, these coefficients are neg- 
ligible insofar as NUq <C 2AE. With this perturb- 
ing Hamiltonian we have not lifted the degeneracy be- 
tween the states with the same occupation of opposite 
wells. If we consider also the perturbing Hamiltonian 



P b 

u j£m u j l tm 



the MS states 
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found in App.lDlarc obtained. These MS states show cou- is not satisfied. We call these coupled states shadows of 
pling to excited states insofar the condition NUq -C 2AE the MS states. 
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FIG. 8: (Color online) Change of the criteria in the high and 
low barrier limit with N . (a) Coupling constant g (dimen- 
sionless) given by criteria (|36[) . (|37[) , (|40[) . and (|43[) for a big 
value of Vq, as a function of N. (b) Barrier height Vb (di- 
mensionless) given by criteria (I34|) , (|35|) . and (|42l) for a small 
value of g, as a function of N. In both cases, the bound- 
aries between different regimes are displaced accordingly for 
increasing N. 
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FIG. 9: (Color online) Eigenj "unctions in the high and low barrier limits, (a) to (c) represent the analytical functions for the 
harmonic oscillator approximation while (f) to (h) are the numerically calculated ones for the Duffing potential, all for Vo = 20. 
The one-dimensional analytical eigenfunctions are shown in (d) and (e) while the numerical ones are shown in (i) and (j). 



